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ABSTRACT 

This paper presents a new analytical treatment of Unsteady Aerodynamics — the 
linear theory covering the subsonic compressible (inviscid) case — drawing on some 
recent work in Operator Theory and Functional Analysis. The specific new results 
are: 

(a) An existence and uniqueness proof for the Laplace transform version of the 
Possio integral equation as well as a new closed form solution approximation thereof. 

(b) A new representation for the time-domain solution of the subsonic compressible 
aerodynamic equations emphasizing in particular the role of the initial conditions. 
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1. Introduction 


This paper presents a new analytical treatment of Unsteady Aerodynamics draw- 
ing on techniques of Operator Theory and Functional Analysis. By Unsteady Aero- 
dynamics we mean here the time-domain solution of the field equations — thereby 
also emphasizing the role of initial conditions generally ignored in the literature be- 
cause of the preoccupation with the oscillatory response. In this paper we consider 
the subsonic (or linearized transonic) compressible case — inviscid of course — in two 
space dimensions. The extension to three space dimensions is in progress. This work 
was initiated as part of the Aeroelastic Stability problem — bending-torsion flutter 
in compressible flow. The incompressible case is treated in [1] for the 2-D strip model 
of Goland (see [3]). 

To clarify terminology, by the “2-D Linear or Subsonic Compressible Case” we 
mean the flow characterized by the partial differential equations (being the linearized 
version of the inviscid TSD, Section 2) for the velocity potential, cp(x, z, t) in the X-Z 
plane, — oo < x < oo and we take 0 < 2 : 


,(1 — M 2 


d 2 4> 
dx 2 


+ < 


d 2 4> d 2 cp 


dz 2 


dt 2 


+ 2 Ma c 


d 2 4> 

dtdx 


( 1 . 1 ) 


The main boundary conditions are 


i) Flow Tangency 

^(z,0 +,t) = w a (x,t), \x\ <b (1.2) 

where w a (x,t ) is the downwash, 2b the strip width, and 

ii) The Kutta-Joukowski conditions (see Section 2). (1.3) 

The initial condtions are 

0 ), ^(x,Z, 0 ). 

The main interest is of course in the acceleration potential 

d(p d(f) 

rp{x,z,t) = —(x,z,t) + Ma 00 —(x,z,t) 
at ox 

(in particular for 2 ' = 0+) and its Laplace transform 

^ /*oo 

tf’(x,z, A) = / e~ xt ip(x,z,t) dt, ReA>0. 

Jo 
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By a “particular” solution of (1.1) we mean a solution of (1.1) satisfying the 
boundary conditions but leaving the initial conditions unspecified — allowed to be 
arbitrary. The first such solution was given by Possio [2] in the form of an integral 
equation for the “oscillatory” case — where the downwash is an oscillation at the 
frequency u: 

w a (x,t ) = w a (x, iu)e lujt , t> 0 
and the velocity potential then is also oscillatory, 

4>(x,z,t ) = z,iu>)e lu}t . 


Several authors (see [3]) have presented various versions of proofs, the clearest perhaps 
being [4] . In this paper we present the Laplace transform version valid for “artibrary 
motion” — for arbitrary downwash functions, where unlike in the oscillatory case, 
close attention has to be paid to the initial conditions; we do this by resurrecting the 
classical time-domain source-doublet integral of Kiissner [5]. New with this paper is 
the existence and uniqueness proof for the Possio integral equation, as well as a new 
closed form analytical solution approximation. Numerical computations using this 
solution for the lift and moment [7] show close agreement with the results obtained 
by series expansions or numerical approximations (e.g., [8]) in the main region of 
interest of the parameters. 

In Section 2 we provide an existence and uniqueness theorem for the time-domain 
solution of (1.1) for arbitrary initial conditions and satisfying the flow-tangency con- 
dition. Whether it satisfies the Kutta-Joukowski conditions depends on the aerody- 
namic initial conditions. Of course it is customary to dismiss solutions which don’t 
as “nonphysical” as in the nonlinear case [9] but the fact remains that in the mathe- 
matical model we have to reckon with their existence. Indeed it is an open question 
to characterize the initial conditions that lead to solutions which satisfy the Kutta- 
Joukowski condition. 

If we set the time derivatives in (1.1) to be zero, we have: 


d 2 cf) d 2 cj) 


+ 


dz 2 


(1.4) 


but if we also retain the boundary conditions, the solution becomes a function of 
time since the downwash is. Again a particular solution can be obtained by classical 
techniques (see [3]) and we shall denote this by z,t). However we should note 

that the initial conditions are not the same for 0 m (®, z,t) and (f)p(x, z,t), the Possio 
solution. If in (1.4) we set 


M = 0; U : a^M ± 0, 
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and U enters via the Kutta-Joukowski condition, we have the “incompressible” case 
which leads to the classical “airfoil” equation [3]. We shall denote this solution 
4>o(x, z, t). It is possible to obtain (p M (x,z,t) from (po(x,z,t ) by a transformation 
of coordinates — the Prandtl-Glauert transformation — but we will not need to 
consider this at all. 

One of our main results is that any solution (p(x, z, t ) of (1.1) can be expressed as 
the superposition 

<p{x,z,t) = <f> M (x,z,t) + 4> R {x,z,t) (1.5) 

where <fi R (x,z,t) is a solution of (1.1) satisfying “homogeneous boundary” conditions, 
but with a forcing term depending on 4>m{') (see Section 5), and also satisfies the 
Kutta-Joukowski if (j>(x,z,t) does. In particular therefore < f> P (x,z,t ), the Possio solu- 
tion, has such a representation. Note that depending on the initial conditions there are 
solutions of (1.1) satisfying the flow-tangency condition but not the Kutta-Joukowski 
conditions — whether “nonphysical” or not. Since the aerodynamic initial conditions 
are never known, this makes “unsteady aerodynamic” calculations (e.g., [10]) difficult 
to verify by experiment. 

In an ascending scale of accuracy we may put 

0o (x,z,t) 

<p M {x, z,t) 

4>p(x,Z,t) 

(all satisfying the Kutta-Joukowski conditions) in the sense that each may be used to 
calculate the lift and moment as a function of U (for 0o('))> or of M, for <Pm(') and 
(pp ( * ) . Note that the initial conditions are different for the latter two, unless we set 

Wa( 0,X) = W a ( 0,X) = 0. 


Organization 

We begin in Section 2 with an abstract (function space) formulation of (1.1) with 
an appropriate (L 2 ) definition of boundary values. We consider first the initial value 
problem for homogeneous boundary conditions leading to the time- domain semigroup 
solution (for appropriate definition of energy) and the resolvent, the Laplace domain 
solution. Unfortunately, some knowledge of the Theory of Semigroups is assumed in 
Section 2. Basically it provides the machinery for going from the Laplace domain 
to the time domain, and in our case both are important. In Section 3 we develop 
the Laplace transform version of the Possio equation, starting with Kiissner’s form of 
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the time domain solution integral. The bulk of the section is devoted to deriving a 
constructive existence (and uniqueness) of solution leading to a closed form solution 
which is accurate to the order 

M 2 log h 

in terms of M or 

|7 2 l°g 7 1 

in terms of 

AM 

7 ~~ U(l — M 2 ) ' 

We show that this solution reduces to the known solution for M = 0. 

In Section 4 we treat (1.4). We obtain a particular time-domain solution, satisfying 
the Kutta-Joukowski condition, for arbitrary downwash function without invoking 
the Prandtl-Glauert coordinate transformation. It is a more rigorous mathematical 
treatment than in [3] with a few new, more general results as well. Finally in Section 5 
we assemble the complete solution — combining boundary value and inital value, 
connecting with the results in Section 2. In particular we show that the Possio 
solution can be written 

Op{i..r.z) = (j) M (t,X,Z ) + <j> R (t,x,z) 

where <f> R ( ■) satisfies (1.1) with homogeneous boundary condition and zero initial con- 
dition but a forcing term depending on <Pm (')■ We also derive an alternate expression 
for the Laplace transform pp(X,x,z) providing an alternate to the Possio equation. 


4 



2. The Field Equations 


The basic Transonic Small Disturbance or TSD equations for the (perturbation) 
velocity potential in compressible flow are, following for example the derivation in [9], 
with <f>(x,y, z,t) denoting the velocity potential: 


d 2 (f) 
d t 2 


+ 2 Ma Q 


d 2 (j) 

dldx 


2 m 9<t> d 2 (f) 2 d 2 (f) 


d 2 4> 


( 2 . 0 ) 


dy 2 00 dz 2 ’ 

— oo < x < oo, — oo < y < oo, 0 < 2 < oo 

where, in the usual notation, M is the Mach number and the speed of sound, and 


£ = ( 1 + 7 )M 2 a 2 00 , 7 being the adiabatic constant. 


In this paper we shall only consider the linear case where we take e = 0; moreover 
we restrict ourselves to the “planar” case where the potential is independent of the 
variable y, so that the partial derivative with respect to y is set equal to zero. Thus 
we have for <f>(x,z,t ): 


.d-M^ + a 


d 2 cp d 2 (p 


dz 2 


dt 2 


+ 2 Ma 0 


d 2 (f) 
dtdx 


( 2 . 1 ) 


The boundary conditions that need to be imposed are: 

i) Flow Tangency: 

(x,0 +,t) = w a (t,x), \x\<b ( 2 . 2 ) 

az 

the right side being the “downwash” function which vanishes outside the interval 
[— 6 , 6 ], 26 being the “strip” width. 

ii) Kutta-Joukowski: 

This condition is stated in terms of the acceleration potential: 


dd)(x,z,t) rT d 

H x,z,t ) = — + U — cf>(x,z,t). 

We require that 

tp(x,0 +,t) = 0 , x = b— and x>b 
(= zero pressure jump at trailing edge). 
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iii) Far-field Conditions: 


0 


(ty = o 

dx ’ dz 

as | a; | — » oo for each nonzero z and z — * oo for each x. See [3], [9] for more on these 
conditions. 

In addition of course considering (2.1) as an “initial value” problem in the time 
coordinate t, we will need to specify initial conditions (at t = 0) as well. 

Abstract Formulation 

The first question to be decided is the choice of the function space in which to 
work. Let R+ denote the half-plane: 

— oo < x < oo, 0 < 2 < oo. 

r 1 3 

Our basic space will be L 2 R 2 + with “weights,” norm (squared) defined by: 

ll/ll 2 = ll/ill 2 + (1 -MV.il/ 2 ll 2 + al\\M\ 2 (2.3) 

where 

fi 

f = h , fie l 2 [rI . 

h 

We shall denote this Hilbert space by 1 H. 


Boundary Values 


Since the domain is not bounded we cannot use the usual “trace” definition for 
boundary values. Here we shall define them in the L 2 sense. We begin with the 
boundary z = 0. We shall say that f(x), —00 < x < 00 , /(•) G L 2 [i? / ] is the 
boundary value on the boundary at z = 0 of a function f(x,z ), |x| < 00 , 0 < z < 00 , 
in L 2 R 2 if 

POO 

/ I f(x)~ f(x,z)\ 2 dx (2.4) 


which is defined a.e. in 0 < z < 00 , goes to zero as z 
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Lemma 


Suppose /(•,•) e L2 R\ and so does y ^ (L2-partial derivative). Then there 
exists h(-) E L2KR 1 ] such that 


/ OO 

I h(x) - f(x,z) | 2 cte 

-OO 


0 as z 


Proof 


We note that 


Tz f{x ’ z) 


dz < 00, a.e. in — 00 < x < 00. 


Let 0 < — > 0 . Then 


/ oo /*oo /• 

|/(®,2n) - /(^,^m)| 2 rfar = / / 

-00 J — 00 J z. 


POO 

P d! f z) dz 

f — 00 

JZm dz 


dx. 


< k 


/ oo /-2: 

I 

-OO J Z n 


9 f{x,z)) 


dz 


dz dx < \z n — z m I 


df 


dz 


Hence f(',z n ) defines a Cauchy sequence in L 2 [R'} which converges to a function, 
denote it h, in L2{R 1 ]- And of course 


/ OO 

\h(x) — f(x, z )\ 2 — * 0 , as z — > 0 . 

-OO 


The limit h(-) is clearly independent of the sequence chosen. Hence h(-) is the bound- 
ary value at z = 0 , and we may use the notation 


h{x) = f(x, 0 +) 


since it would be consistent with the pointwise limit, should there be one. 

The situation is different on the boundary x = +00, x = —00 in that if, as we 
shall, we define the boundary value to be /(+00, z), if 

\f(x, z) — /(+00, z )\ 2 dz — » 0 as noo 

(and similarly /(— 00, z) as the boundary value at x = —00, if 

\f(x, z) — /(— 00, z )\ 2 dz 0 as x —>■ —00 ) 
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in that 


/x(v) ^ L 2 [/hi 

does not assure this. Although we do have from 

fi ,9 


|/(L,z)| 2 = 2R e^ — / (x, t) f (x, z) dx - |/(0, 2 )| 2 


and f & L 2 


i? 2 


that 


0 

0 


as L — > oo 
L — — oo. 


as 


f(L,z) 

f(L,z) 

Similarly we define the boundary value / (x, oo) by 

\f(x,oc) - f{x,z)\ 2 dx 0 


a.e. m z . 


/: 


as 2 


oo 


(2,5) 


assuming of course the limit exists. 

We note that (2.1) is an initial-value plus boundary- value problem, albeit linear. 
The technique for solving such problems using the Theory of Semigroups of Operators 
is outlined in [7], Thus we can construct the solution as part due to the initial value 
problem setting the boundary value to be zero and part due to the boundary problem 
without regard to the initial value. The former leads to a Cauchy problem and that 
is what we shall deal with first. 


The Initial Value Problem: Semigroup Solution 

We introduce the operator A with domain and range in 7i by: 

and 4A G A 




fl 

V(A) = 

f = 

h 



h 


dx 


dz 


a£ and §7 e L 2 


R\ 


Rl 



with the “homogeneous” boundary conditions: 


\f 3 (x,e)\ 2 dx 


0 as £ 


0 


and the “regularity” conditions: 

/“OO 

\fi(L, z)\ 2 dz —> 0 as \L\ 


n 


oo, i = 1,2,3, a.e. in x 



/ oo 

\fi(x, L)\ 2 dx — » 0 as L — ► oo, a.e. in x 

-OO 


A/ 


-2t/£ 


_9_ 

dx 


d_ 

dz 


0 


0 


fl 


h • 


/s 


( 2 . 6 ) 


Note that the Kutta-Joukowski condition is omitted. Thus defined, A has clearly a 
dense domain. Moreover: 

Lemma 

A is dissipative (on its domain). 

Proof 

With [ , ] denoting inner products in Li. as well as L 2 [/?/], we have: 


WJ\ = 


-2 Udf 


1 , /i i/2\ 2 df2 2 U J 3 n 

- + (1 - M )a tx — + a iX> — , fi 


dx 


<9/3 


dx 


+ (1 - M 2 )a 


2\„2 

'OO 


df 


and 

2 Re [AfJ] = -2 U 


dfi 


i , fi 


dx 

+ (1 - M 2 )a ; 

dfi 
dz 


2 \ „2 

’oo 


+ 

d_h 

dx 


f 


i ? 


5a: 


5A 

5a: 


h 


dz 


+ « 


dfi ' 


fl 


+ 


5 /3 


+ 


/s 


/l, 


dfi 

dz 


dh 

dx 


+ 


dfi 

dx 


, /2 


dh 

dz 


fi 


f ^ 

.72 j q 
OX 


' ^ 
J 1 ? ^ 
oz 
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Now 


r fL fdf!(x, Z ) — , t( x df 1 (x, z) 

Jo dz J. L {—a ~ ~ h ^ z) + h ^ z) lh~ 

J Q (\M L > Z )\ 2 “ \fi(~ L i z )\ 2 ) dz 


dh ,1 \r d fl ] \d fl 1 r dh\ 

"TJ , Jl + / 2 , TT - + “T - , J 2 + Jl , 

OX ox ox ox 


= lim / 

L->oo JO 


f 2 (x,z) fi(x,z) + fi(x,z) f 2 (x,z) dz 


by virtue of the boundary conditions. Next 


df-i d f r 00 -| 

wjr , / 3 + /i, wjr = lim J [f 1 (x,z)f 3 (x,zj\dx 

J L J L— > oo 00 

£->0 + 


Similarly 


§M + f/,.^1 = o 


where we have used the vanishing boundary conditions at x = ±oo and z = oo. As 
for z = 0 we note that 

/ oo _ / r oo / /‘oo 

fi(x,z)f 3 (x,z) dz < J \fi{x,z )\ 2 dx J / |/ 3 (x,z)| 2 dz 

-OO V J — OO V J — 00 

and as z — > 0+ since ^ G L 2 [i? 2 ] 


|/i(2:, z ) | 2 dx 


|/i(2:, 0+) | 2 dx < 00 
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and 


/ 1/3(2;, z )\ 2 dx — » 0 as z —> 0 + . 

J — OO 

Hence it follows that A is dissipative and has a dense domain, and A can be closed 
with the closure defined by 

A = (A*)*. 

But using the usual arguments involving distributional derivatives, it is readily seen 
that: 

A* = —A on the domain of A 


A* = -A; (A)* = -A. 

Hence A is dissipative and hence is the generator of a contraction Co-semigroup which 
we shall denote by S(t), t > 0 . We hasten to remark that 

(Sit))- 1 = s(ty 

and hence S(t) becomes a group by defining 

S(~t) = sit )- 1 


and of course 

||<f?(£)|| = 1 for each t, —00 < t < 00. 

What is important for us is to note that the point spectrum of A is empty. (A 
has no eigenvalues.) The resolvent set of A is the set of those complex numbers A 
such that the nonhomogeneous equation: 

Xg ~ Acj = / 


has a unique solution for every / in H, and we use the notation 

g = r(\, A) f. 

From the dissipativity of A, it follows that the resolvent set of A includes all A such 
that Re A > 0 , and we actually have the construction 

R(\A)f = J°° e~ xt S(t)f dt, Re A > 0. 

In fact it can be recognized as the Laplace transform version of ( 2 . 1 ). Let 



fl 


9 1 

/ = 

h 

, 9 = 

92 


fs 


93 
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Then (assuming that / is such that g is actually in 'D(A ) ), we have 


9i = 

X 4>{x, z , A) 


d ~ 

92 = 

a-J (x ' z ' x) 


d ~ 

9s = 

— 0(x,z,A) 


where <f>(x,z, A) is the solution of 


\ 2 <fi(\,x,z) + 2Ma 00 X 


d(f>(\, x , z) 


a 


dx 

2 d 2 (p{X,x,z) 
dz 2 


al(l~M 2 ) 


d 2 (f>( A, x, z) 
dx 2 


= fi(X,X,z), — oo < X < OO, 0 < 2 < oo. 


Resolvent 


To evaluate the resolvent, it is convenient to work with Fourier transforms. Thus 


for /(•, •) in L 2 




we define (in the usual sense): 

^ r oo /*oo 

f(iu> i , Lj 2 ) = / / e ~ luJlX ~ luJ2Z /(x, z) dx dz, —00 < , o ; 2 < 00. ( 2 . 7 ) 

J — oo JO 

Also since for cr > 0 

/ oo / /*oo \ 2 />oo 1 /*oo 

( / e _CTj: |/(x, £)| ete ) dx< — e~ az \f(x, z)\ dz dx 

00 \ J 0 / J— 00 O’ JO 

^ 7*00 /*00 

< — / / |/(x, 2:)| 2 d.2 dx < 00, 

(T J— 00 J 0 

we can also define the Fourier- Laplace transform 

/ oo . /-oo 

e lWlX dx / e A iz f(x,z)dz , Re/i>0. (2.8) 

-00 J 0 

In particular we see that for / in 12(A), 


f 


fi 

h 

h 


/i > 0, 
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Jirr^ ///] (cji , //) = / e ia;i */i(:E,0+) ah 


which we can also express as: 


Let 

Then noting that a.e. in x : 
dfi(x,z) 


l '■ 

— / fi(ioJi , hn 2 ) duj 2 

Z7T J — 00 

9 = Af. 


-IUJ2Z 


dz = -/i(x,0+) + iuj 2 e tW2Z f 2 {x,z) dz 
az Jo 


we have: 


9 = H 


fi 

h 

h 


0 

0 


/ oo 

e -iuixf^ x ^ o_|_) dx 

-OO 


where H is the multiplier matrix: 

—WiuJi a 2 ^! — M 2 )(iuJi) a^iu 2 


JT 




2^2 


0 


0 


0 


0 


Note that we can use (2.11) to define A as a “multiplier” operator. Moreover 

A f ~ Af = g 


becomes 


9 = (A Is - H)f + 


0 

0 

e ~ lu}lX fi(x,0+) dx 


and hence 


/ = R(\,A)g, ReA^O, 


(2.9) 

( 2 . 10 ) 

( 2 . 11 ) 

( 2 . 12 ) 
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becomes 


- h r 1 


(XIs - H )~ l g. (2.13) 


/ OO ^ 

h(iuJi , zcu 2 ) dcu 2 

-OO 


Noting that 



A 2 

XciiuJi 

Ac 2 ic^ 2 

= A) 

Aicui 

A 2 H - 2A U iuj\ -|- c 2 cl?2 

c 2 o; 1 o; 2 


— Xiu-2 

c 2 o; 1 o; 2 

A 2 -f- 2A U iuj\ CiCcJ 2 


where 

Cl = (1 -M 2 )a^, c 2 = 

rf(A) = A ^A 2 “I - 2 £/ ic^iA -|- cic<j 2 -|- c^cj^ 

7 ^ 0 for A ^ zcj, a; real , 

we have, writing /r in place of iu 2 : 

A 

/ = A , 

A 

- A 2 <)i + Ac] (iuJi)g2 + A/m 2 A 3 — Aft) 

f 1 V7v2 i o771 \ i 1 32 oT j Re /i ^ 0 

A (A^ H - 2U VjJ\\ H - Cic^i — c 2 ^x ) 

where fib is a function of hni defined by 

/ib(hni) = lim^ /i/i (tcui , g) 


(2.13a) 


A = (a^!^ + (A 2 + ‘IXlJiuJi - c 2 g 2 )g 2 - c 2 g(iuJi ) (</ 3 - As)) , Re/i > 0, 

(2.13b) 
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h 


d( A) 


(^Xfigi + c 2 {iu} 1 )ixg 2 + (A 2 |- 2\Uiu\ + c±u^) (jj s 



We should also note the relationships: 

fi{x,z) = A [ f 2 (y, z) dy 

J — oo 

and when g 2 (y, 0+) is defined: 

fi{x,0+) = A / f 2 (y, 0+) dy 

J — OO 


[ g 2 (y,z)dy 

J — oo 



#2(^,0+) dy. 


Re //, > 0. 
(2.13c) 


We know of course from the general theory that i?(A, A ) is defined for A not pure 
imaginary. But (2.13a)-(2.13c) may well continue to be defined for A imaginary as 
well for some g since the range of 

(iul - A ) 

is dense in 7i. Indeed if 

d( A) = 0, 

then 


A 


-2 Uicj! ± s J-4ujU 2 - 4( Cl uf + c 2 u$) 

2 

i —Uui ± \/u) 2 U 2 + Cicof + c 2 u>2 ■ 


We need only require that the numerators in (2.13a)-(2.13c) are also zero for this 
value of A, so that /,; defined therein with iuj 2 in place of /i are square-integrable. 


15 



Spectrum of A 

The resolvent of A is clearly not compact being characterizable as a “multiplier” 
operator. Hence we cannot expect that A has eigenvalues. In fact the point spectrum 
of A is empty, and it has a pure continuous spectrum: Re A < 0. Again this follows 
from the “multiplier” characterization. Indeed 

Af = A / 


yields from (2.13a) that 


fifai , iu 2 ) 


All>2 C 2 fib 


(A 2 T 2zAiiqA T Ci^i T ^ 2 ^ 2 ) 
and integrating with respect to u >2 yields: 

^ POO ^ 

fib = / fi(ioJi , ito 2 ) du 2 = 0 . 

J — OO 

Hence f\ is zero also, and so is /, in turn. 


Velocity Potential 


Given the semigroup solution 


F(t) = S{t)F{ 0), F(0)eV(A) 


where F(t) denotes 

fi(t,x,z) 

h(f,x,z) 

fs(t,x,z) 

we can relate the solution to (2.1) in the following way. Define the velocity potential 
(f>{t,x,z) by: ^ 

(f>(t,x,z) = 0(0 ,x,z) + / fi(a,x, z,) da (2.14) 

Jo 

where 0(0, x, z) is required to be such that 


<90 <90 

dx ’ i)z 


Lo 


Rl 
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Then we can calculate that 


d 2 0 w d 2 4> 

dt 2 dtdx 


(1 - M 2 )a 


2 

oo 


d_h 

dx 


+ o, 


dfs 

dz 


Hence 


We define 


so that 


Similarly 


Hence 


We define 


d 2 4> 9 2 0 dj\ df 2 

dtdx dxdt dx dt 


d 


dcf) 


d 2 (p 

dtdz 


dt 

| dx 

90(0. 

,X,z) 

dx 


dcf) 


dx 

1 

d 2 cf) 


dzdt 

9 

90 

dt 

dz 

90(0. 

,x,z) 


dz 



= o. 


= f2(0,X,z) 


— h ■ 


dfi = dfs 
dz dt 


~ fs 


= 0. 


: fs(0,x,z) 


and hence 

d(j) 

dz 

Hence it follows that 0(t, x, z) satisfies 


h ■ 


with 


, ,,, 

dt 2 dtdx 


(i 


M 2 )a: 


d 2 cj) 

dx 2 


T 


d 2 (j) 
dz 2 


0(0, x , z) given 


d 

dt 


<f>(t,X,Z ) 


t = 0 


given 


fi(0,x,z) 


(2.14a) 


(2.15) 

(2.16) 

(2.17) 
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t=0 


/ 2 < 0 , x, z) 


(2.18) 


dcj) 

dx 


(t : x,z,) 


where 


dcj) 

dz 


(t,x,z,) 


t = o 


fs(0,x,z) 


F(0 7 x,z) 


fi(Q,x,z) 


f2(0,X,Z ) 

fa(0,x,z) 


<E V(A). 


(2.19) 


Also we see that the solution is unique, subject to the initial conditions (2.16)— (2.19) 
and the boundary conditions imposed in describing the domain of A. Note also that 
from (2.14a) we have 


Hence 


dfz d dcj) d dtp dfi 

dt dt dx dx dt dx 


fl(t,X,Z 



0 f 2 {t,y,z) 

dt 


dy 


( 2 . 20 ) 


Kutta-Joukowski Condition 


So far we have not addressed the Kutta-Joukowski condition. In view of (2.14) 
where we may think of fi(t,x,z) as the time-derivative of the velocity potential, and 
f 2 (t,x,z) as the derivative 

dcj)(t, x , z) 
dx 

we may state the Kutta-Joukowski condition as 


_/j(f, x, 0+) + Uf 2 (t,x, 0+) = 0 at x = 6— 

= 0 at x > b. 


( 2 . 21 ) 
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We may then include it as a restriction of the domain of A: 


D(A 0 ) C D(A) 


D(A 0 ) = [f | / e T>(A) and / satisfies (2.21)] 

Aof = Af. 

The restriction Aq thus defined still has a dense domain and is of course dissipative 
thereon but the closure of Ao is readily seen to be A. 

We can also give a “dynamic” Kutta-Joukowski condition in the following way. 
The acceleration potential defined by 


ip(t,x,z) = 


d<j>(t,x,z) dcj)(t,x,z ) 

dt dx 


can by virtue of (2.20) be expressed in the form 


'ip(t,x,z) = I df2 % V ' ^ dy + Uf 2 {t,x,z). 


( 2 . 22 ) 


The Kutta-Joukowski condition can then be expressed entirely in terms of f 2 (t,x,z) 
as a “dynamic” condition since it involves the time-derivative. Thus ^>(t, x , z) defined 
by (2.22) is required to satisfy 


0+) = 0 at x = — b and for x > b. 


We shall refer to this as the “dynamic” Kutta-Joukowski condition. 
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3. Solving the Boundary Value Problem: Possio’s Integral Equation, 
Laplace Transform Version 


In this section we obtain one solution (j){t,x,z ) to the boundary- value problem: 


d 2 (h d 2 (b 

— + 2 U 

d t 2 dtdx 


al(l-M 2 ) 


d 2 0 

dx 2 


+ O'n 


d 2 (j) 

dz 2 


(3.1) 


satisfying the flow tangency condition, the Kutta-Joukowski condition and the far 
field conditions, but without regard to any initial conditions. 

We begin with the construction for a particular solution of (3.1) due to Kiissner 
[5, p.7 ] (see [18], for a direct proof that (3.1) is satisfied), using doublets at the 
source: 


</>(x,V,z,t) 


| rb roo /»£— £ 

— / dt; / dr] dx' 
4 l 7 T J—b J — oo J—oo 


d_ 

dz 


A (f, r/, t + 


x' _ x-i _ \/ x' 2 +(l- Af 2 )((^- 97 ) 2 + 2 2 ) 
( 7 ( 1 -A/ 2 ) U aooil-M' 2 ) 


\Jx' 2 + (1 - M 2 )((y - T]) 2 + z 2 ) 


where A(£,r],t) is the doublet intensity on the airfoil which is to be determined from 
the boundary conditions. Note that 


4>( x i Vi —z, t) = -0(x,y,z,t). 


And further, 

<p(x,y, 0+,0+) = 0. 

Since cp(x,y, z,t) does not, in our case, depend on y , we have: 

^ rb roo /*#— £ 

(f>( x ,z,t) - ~r / / dr i dx> 

47 T J — b J — oo J — oo 

i 

U ' ( 7(1 — A/ 2 ) 


d_ 

dz 


Alt f , X 1 _ aA ' 2 +(1 -M 2 )((>? 2 +z 2 ) 

^ ' IJ ' TKl-M 2 ) aoc(l-Af 2 ) 


\Jx' 2 + (1 — M 2 )(r] 2 + Z 2 ) 
where A(^rj,t) does not depend on 7] either, so that we may denote it A(£,t). 


(3.2) 


Initial Conditions 


We note that (3.2) defines the initial ( t = 0) values of the velocity held. It is easy 
to verify from 


A(£,t) =0, t < 0 
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that 


(f)(x,z, 0) = 0 
(j)(x,z, 0) = 0 


from the fact that the set 


x' < x — £ 




+ 


X 


\Jx' 2 + (1 —M 2 ){rj 2 + z 2 ) 


U ' U (1— M 2 ) 
is actually empty, for M < 1. 


a^l-M 2 ) 


> 0 


Since 


d 


W a (x,t) = —(j)(x,z,t) 


2=0 


we have from (3.2) that 

^ pb poo 

w a (x,t ) = — I I dr) dx' 

47 T J — b J — oo J — 


6> 2 


^11 _|_ x' _ ~ y / x' 2 + (l-M 2 )(r]' 2 +z' 2 ) 


U 1 U(l^\l 2 ) 


Qoo((l -AT 2 )) 


+ (1 —M 2 )(r) 2 + z 2 ) 


(3.2a) 
1 2=0 


Unfortunately this is as far as we can go in the time domain. We therefore invoke 
either Laplace or Fourier transforms. Traditionally the Fourier transform is preferred 
because one may think of an “oscillating” doublet corresponding to an oscillatory 
downwash input. The first such derivation is due to Possio [2] leading to the Possio 
integral equation, the most lucid derivation of which may be found in [4] . As we show 
later, the Laplace transform can formally be obtained from the “oscillatory” version, 
and also vice versa. Here however we shall actually derive the Laplace transform 
version of the Possio equation, emphasizing the role of initial conditions. It will turn 
out that this is actually more useful than the time-domain solution since generally 
the primary interest is in the stability of an aeroelastic system [1], A Laplace trans- 
form version is given in [11, p.3] but it involves divergent integrals which must be 
interpreted carefully. In fact the author derives a formula ([11, p.4, eq.31), for the 
initial conditions therefrom which is invalid in the present (2-D) case. 


21 



We note that 

Afar) = 0, r < 0. 

Let 

x-£ x' \Jx' 2 + (1 - M 2 )(rf + z 2 ) 

~U U{ 1 - M 2 ) + Ooo(l - M 2 ) 

Then since we must have that 

t — a > 0 for all t > 0 


we need only consider 

a < 0. 

Hence defining the Laplace transforms: 

-A t 


A(£, A) = / e M A(£,t) dt] <j>{x,z, A) = j e xt 4>(x,z,t) dt, ReA>0 


we have, taking Laplace transforms in (3.2): 
1 r 6 /**-£ 


* 1 /■0 

4>( x ,z,\) = — I f d t da: 


exp A 


:r: 


s - 

H(1 - M 2 ) C/ 


i(e,A) 


( \ \J x' 2 + (1 -M 2 )(rj 2 +z 2 ) \ 

a f 00 exp V " ^ acc,(l-M 2 ) j 

Oz J-oo (1 —M 2 )(r] 2 + 2 2 ) 


dry. 


Making the change of variable, 

(1 — M 2 )r/ 2 


- 1 


we have: 


C > 1, c = 


2 a 2 + {l-M 2 )rj 2 


a 2 = x' 2 + (1 - M 2 )z 2 


\jx’ 2 + (1— M 2 )(rj 2 + 2 2 ) 


cr 


dr) 


°L i dC u 2 = 

1 — M 2 r/ ^ ’ 1 (1 — M 2 ) 
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( \ ^ x ’ 2 + ( 1 -M 2 )(j) 2 +Z 2 ) \ 

GXP ^ A aoo (l^f 2 ) J 

x ' 2 + (1— M 2 )(rf + z 2 ) 


dr) 


r™ / -A a( \ J_ C « 2 

Ji exp l aoo {i-M*)) ' C« ‘ n ' (i -M2) 




0 r / - A < 

z / exp 

7 l \ ’-"OO 


1 Cv 7 ! - M 2 


a 


a^(i-M 2 ); c« avc^n: (i-aa 2 ) 

/ -AaC \ 1 


<*c 


Vl - M 2 Ji eXP Va 00 (l-Af 2 ) y / 


d< 


Vl - M 2 


tfol 


X\Jx ' 2 + (1— M 2 )z 2 
aoo (l-M 2 ) 


Hence 

4 >(x,z,X) = j_ b di • i(£,A) 


L 


X-S ( 1 


A / ad 


vTUh 


exp 


C/Vl- Af 2 


0 , (X y / X ' 2 + (1-M 2 )Z 2 

^ Ao ( a^j 


fa -0 




(3.3) 


which is our basic formula in what follows. Note that: 


9 ^ ^\J x ' 2 + ( 1 —M 2 )z 2 

dz \ a co (l-M 2 ) 


-IV 


X \Jx ' 2 + (1 —M 2 )z 2 \ A 


(1 -M 2 )z 


a 00 (l-M 2 ) J a 00 (l-M 2 ) ^2 + (1 -M 2 )z 2 


X ( A \Jx ' 2 + (1-M 2 )z 2 

O'oo ^(l-M 2 ) 


We invoke next the boundary condition 


x ' 2 + (1 — M 2 )z 2 


(3.3a) 


d 

W a (x,t ) = — <(>(x, 0 +,t)- 
cJz 
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Denoting the Laplace transform 


and noting that 


/» OO 

/ e~ xt w a (x,t.) dt. by iv a (x,X) 
Jo 


roo Q Q 

j e~ xt — (f>(x,0+,t) dt = — 4>(x, 0+, A), 


we have from (3.3) that 

1 f'°° - 1 r x ~Z ( A / v 

= 2? exp £/(tVm5 


- (s-f) 


t) 2 ^ /A^ + (1_ M 2)^2\ 

a oc (l-M2)) j 


or, alternately, 


- , „ I* .4«. A) rfi; 

V ’ ' +J- b 2tt VI - M 2 


/•*-* /-A , A 

L exp (— <;r - f _ ^ + F 


A / M 2 
U V V 1 - A-/ 2 


Ko ( a , x (i-j/^) \ /s ' 2+(1 ' M2)22 ) d!/ 


lim — [ b A(€, A) ^ 

2^0+ 27T J-b Vl — M 2 


(- X(x-0 

o exp — tr- 


■ exp (f (IFSP)) FF A °( a x (l'-JI/ 2 ) vV + (l-M 2 )z 2 )] 


where we note that 


3z 2 Ao(,) <&(1-M 2 ) A °(o x (l-.V 2 ) \A 2 + (1 V 2 )z 2 


- ^ $ A »(zm^ ^ 2 + (1 - M2)z2 
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where we have used the identity [12, p. 79] : 

-K\(x) 


K[(x) 


x 


- Kq( x ) 


K( x ) = -Ki(x). 


Relation to Pressure Distribution 


The acceleration potential is defined by 


dd> deb 

= - + U-. 


Denoting the Laplace transform 


we have 


since 


e M xp(x, z,t) (it by tf(x,z, A) 


r) - 

ip(x, z, A) = \(f>(x,z,\) + U — 4>(x,z, A) 

ox 


<t>{x,z, 0+) = 0. 

Using (3.3) for (f)(x,z, A) we obtain 
U 1 rb 


^(x,z,X) 


2 tt Vl-M 2 J- 


//(?, A) 


exp 


~M x ~0 + A ( x - 0 


U U (1— M 2 ) 


K ° ( TTTWTm 


d 


dz \U(l—M 2 ) 


U 1 


2 tt Vl-M 2 J-b 


/1(U A) exp 


-A M* 

^(.t-0 2 + {l~M 2 )z 2 \ u 


U (1 -M 2 ) / 

A ^/(ar - C) 2 + (1— M 2 )z 2 


We see that for any e > 0, 

U 1 

lim DnLr, A) = lim — , 

cD 0V ’ ’ ' *-o 2 tt v / 1—M 2 


-A 

TT 


(1-M 2 ) 
'A(x-C) M 2 


L a </ 4 A( t’ x> exp ( : u ( i - jw 2 ); 

MA ^/(a; - £) 2 + (1 -M 2 )z 2 


(x-£) 2 + (1-M 2 )z 2 


U 


(1 —M 2 ) 
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for any e > 0. Further for e and z small enough, 


can be replaced by 


MX a J(x- 0 2 + (1 - M 2 )z 2 

~U (1— M 2 ) 


U (1— M 2 ) 

AMJ(x-0 2 + (1-M 2 )2 2 


Hence for £ and 2 small enough: 


= / 
1 —M 2 J 


(1 -M 2 )z 


-S a 2 + (1 -M 2 )z 2 


da A(x, A) 


u r- d 


<-»£/ 


da 


zx/l^M 2 


Hence it follows that 


Vim^(x,z,X) = (~l) — A(x,X). 


The pressure distribution defined by 


AP{x, A) = ( p) lim 2 , A) - ^(x, - z , A) 
= (~P) lim U’(x, z, A) - $(x, - 2 , A) 

2 J 

= C/i(x, A )p. 


Hence 


A(x , A) 


AP(x,\) 


In particular, by the Kutta-Joukowski condition we must require A(x, A) is continuous 
at x = b— and 

A(b~, A) = 0. 
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Possio’s Integral Equation 


To derive the Laplace transform version of the Possio integral, we begin with (3.4) 
where 


r x-£ 


exp 


A y \ d 2 


U(l—M 2 ) J dy 2 ya 0 (l— M 2 ) 

x-i 

exp 


K 0 \ 

Ay 


A 


sjy 2 + (1-M 2 )zA dy 


Al<o\ 


X 


U{l-M 2 ) dy \a 0 (l— M 2 ) 


\Jy 2 + (1—M 2 )z 2 


X r x ~Z 


l 


exp 


A y 


U (1— A/ 2 ) J-oo 1 \U(1-M 2 ) 


( X ( x -®\ rl 

~ exp yu (T^m 2 ) ) ,Ai1 


A 


T-(l -M 2 ) 

A (x - 0 


(x - £) 2 + (1-M 2 )z 2 


M 1-M 2 ) ^ (,r -^)) 2 ; (\ -M 2 )z 2 


X 


x-£ 


U{l-M 2 ) 


/ % \ r ' / A 

eX P ( TTt, I A o( 


17(1— M 2 ) I \a 0 (l— M 2 ) 


^-0 2 + (i-M 2 )^ 2 


A 2 /•(*-« / Ay 

exp 


U 2 (l—M 2 ) 2 


U{l-M 2 ) / 
X 




which as 7 — >■ 0 


/ A(.r - 0 \ ^ 
ex P ( — — 7777 I R i 


X 


U(l-M 2 ) \ aoo (l—M 2 ) 


\x -f| 


A — — 0 


a QO (l-M 2 ) |a:-£| 


A ( X(x-0\ T , i 

ex P ( — — 777“ I An | 


A 


17(1— M 2 ) ^ V LA(1— A/ 2 ) / \ a QO (l— M 2 ) 


|x - f| 


A 2 


u 2 (i—M 2 y 


■x-£ 


exp 


Ay 


/L 


A 


U{1-M 2 )J u Va QO (l -M 2 ) 


|y| dy 
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Hence (3.4) can be expressed: 


J_ r» A{£, A) 
27T J-b y/l-M 2 


exp 


-Apr -O ' 
U 


(— 1)(1 -M 2 ) ( ex p4 !', f — ( K i 


A 


k-CI 


A 


— (x 


U(l-M 2 ) V V^a-M 2 ) 1 s 7 a^l-M 2 ) |.t- 

A / A 


AT 0 | 


® - f| 


+ 


A 2 


U (1 — M 2 ) u \ aoo (l-M 2 ) 1 

Ay \ r , / A|y| 


U 2 (l—M 2 ) 2 


L r p 


U(l—M 2 ) 


I< a 


a^l-M 2 ) 


dy 


+ 


A 2 


a oo(l — A4" 2 ) J- 


( x ~t) ( A y ( A|y| 

ex P — — A o| 


U(l—M 2 ) u \aoo(l -M 2 ) 


dy 


i * i(e, a) 


27r J-b Vl -M 2 



A 


aoo(l-A'/ 2 ) 


' X(.r-C)M-\ 

, exp 77(1^ J 

A ± 2 _« + A x„( — ix 


A 2 /-M/ {-Hx-0 

exp - — - + 


aoo \x-£\ U "\ aoo (l -M 2 ) 
Ay \W / A|y| \ 


U 2 J-c 


u 


U(l—M 2 ) 


Kn 


««,(! -M 2 ) 


dy 



Mt, A) 

Vl-M 2 


^6 



A (x-QM 2 \ 
U(l—M 2 ) ) 


A i 


AM 

U(l—M 2 ) 



a M[x-ei 

A (rr-0 


A 2 

H 2 


rO-£) 


exp 


Ay 


U{l—M 2 ) 


A (z ~g) ' 
U 



AM|y| \ 

u (i— m 2 ) / 


i^o / AMlx-d \\ 

(1— M 2 ) y[/(l— M 2 ) yy y 



We can simplify this a bit further by noting that [13, p.708]: 


/l( exp w£=m) Ko{ 

Hence hnally we have: 


Alyl \ , u rr^ to, + 


W a (x,\) = J b G (\ x-£)A(£, A) d£, \x\<b, 


where the kernel G'(A, •) is defined by: 
1 A 


G(\,y) 


Vl-M 2 U 


M\y\ /A M . A /A M 2 ^ 

^ll 77 77 7777 |2/| exp - „ ~ V 


y \ U (1 — M 2 ) 
A M 


U (1—M 2 ) 


+ K 0 Ur 


U (1— A/ 2 ) 


I A /A Af 2 
M ex P 77 77 7777 2/ 




7, a|«7|m\ 

/ mv7j exp 


U (1—M 2 ) 

— Ay An 


+ 


(3.8) 


U U(l-M 2 ) 


r— 77, / 1 + Vl^M 2 \ /-Ay' 

log — — exp 


M 


U 


da 


(3.9) 


The integral in (3.8) has again to be interpreted in the Cauchy sense. 


Remark 

For M = 0, (3.8) reduces to 


T’ a (x, A) 


T J* A(t,\)G(\,x-S) dZ 


where 


G(A,y) 


1 

y 


A 2 

H 2 



log 


y-y 1 

|y| 


da 


1 

y 


A 

t/ 


Chi 


Ay 

c/ 


Shi 


Ay 

u 


For the definition of Chi (•), Shi (•) see [13]. This is a known result for the “oscillatory 
case”: A iuj ; see [3,4,5]. Note that unlike the “oscillatory case” we do not use 
divergent integrals and “retain only finite parts” as in [4] . 
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Remark 

We do have that 


Xw a (x, A) 


d , x 
U — w a (x,\) 

dipo(x, 0, A) 
dz 

U 


1 


lim 

z ->° 2n sJl-M 2 


A) 


exp 


A(x — £) M 2 


£/ (1-M 2 ) 


a 2 


K 0 


AM 


dz 2 \U(1—M 2 ) 


(x-0 2 + (l~M 2 )z 2 


which we could use in place of (3.8) since the left side is given, but we cannot take 
the limit inside (in the integrand) because it leads to a kernel with a singularity of 

1 

(;r - f ) 2 

and the corresponding Cauchy integral is not defined! This is even more evident for 
M = 0. Solving for w a (x, A) essentially “smooths” this to make the singularity of 
order one: that is, 

1 



x - f ' 


Possio’s Equation: Oscillatory Case 

We can obtain the “oscillatory” case from (3.8) as follows. Let 


Then 


w a (x,t ) = w a (x,iu>o)e lulot , t > 0, <n 0 real ^ 0. 

w a (x,iu 0 ) 


w a (x, A) 


A — iojQ 


Hence (3.8) yields 


w a (x,iu 0 ) 1 f b c\ AC 

— : = — A (£,X)G(X,x -f ) dC 

A — tU 0 Z7T J-b 


Hence the right side has a simple pole at A = icu 0 . We note that G(X , •) has no poles 
on the imaginary axis, if we omit zero. Hence 


i(£,A) 


^o) 

A — iujQ 
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Substituting and letting A — > iu 0 , we obtain 

1 f b - 

w a (x,iu 0 ) = — A(^ } iu 0 ) G(iu 0 ,x - £) d£. (3.10) 

Z7T J-b 

This is the Possio equation. See [3,4] where Hankel functions are used in place of 
the Bessel K functions. Conversely given (3.10), we may formally replace iuj by A to 
obtain (3.8), provided the kernel is analytic in the right half-plane. In other words, 
we may replace iui by A in the oscillatory Possio equation, and we could get (3.8) in 
this manner from the oscillatory versions in [3,4] using the Bessel K functions. 


Existence and Uniqueness of Solution 

Uniqueness 


We begin with Uniqueness. Suppose then that for some nonzero A, Re A > 0, (3.8) 
has two solutions in L ^~ £ , for some e, with the property that each vanishes at 1— . 
Then the difference, denoted A(- , A) will satisfy 


0 = 




hi < b. 


(3.11) 


Going back then to (3.3), <f>(x, z , A) will satisfy the Laplace transform of (2.1), 


(3.12) 


Moreoever the 


satisfies the boundary condition 
dtp 


dz 


(. x , 0+, A) = 0, hi < b. 


Hence as we have seen 


F 


satisfies 


A F 


\<p 

d<j> 

dx 

dtp 

dz 


AF. 


but the point spectrum of A is empty. Hence F = 0, and </>(•) = 0, proving uniqueness 
of solution to (3.8). In other words the uniqueness is a consequence of the fact that 
A has no eigenvalues. 
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Existence 


The question of existence (or solution to (3.8) for each A, Re A > 0) is not an 
idle one — the existence usually comes with some means of constructing the solution. 
This is indeed the case here, as we shall show, drawing on the work of Tricomi [14] 
and Sohngen [15] on the airfoil equation. 

We begin with some Lemmas. 


Lemma 1 

Let / G L p [— b, b\ ( L p for short in what follows). Then the finite Hilbert transform 
H, defined by: 

Hf = g 

g(x) = \f 1 — M ' 2 — I a.e. |rc | < b (3.13) 

where the integral is defined in the Cauchy sense, defines a bounded linear operator 
on L p into L v for any p > 1. 

Proof 

Titchmarsh [16, p. 132] . 

Lemma 2 


Let g(-) £ L°°[— 6, b] ( L°° hereinafter). Then 

Rg = f 

2 \ i 


f(z) = 


X 


irVl -M 2 J V b + x 


b b + £ g(0 


b~£ 


a.e. |rr| < b 


(3.14) 


defines a linear bounded transformation on L°° into L3- 


Proof 

Sohngen [14] , Tricomi [15, p.175 et, seq.]. 
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Lemma 3 


Suppose g 6 L°°. Then the integral equation: 

/// g 

is satisfied by any /(•) of the form: 

f(x) = ( Rg)(x ) + ^ 

where c is an arbitrary constant. Furthermore if Rg is bounded at x 
is the only solution with that property. 

Remark 

(3.15) enables us to construct solutions of (3.12) which do not satisfy the Kutta- 
Joukowski condition. 

Proof 

Sohngen [14] , Tricomi [15]. 

Corollary 

As an application of Lemma 3, consider the example (setting 6=1): 

1 p df 1 

g(x) = toLjTT) = 2^ log ^ 1 + ^ ~~ log(l-z)), M < 1- (3.16) 

It follows that 

( Rg)(x ) = 1, — 1 < x < 1. 

Let 

g+ (x) = log(l + ar) 

9-( x ) = -log(l-x). 

Then straightforward analysis shows that 

(. Rg + )(x ) — >• 0 as x — > 1. 

Hence it follows that 

(Rg )(x) — ► 1 as x 1. 
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Lemma f 


Suppose (](■) satisfies a Lipschitz condition for some a, +| < a < 1 


\g(t) — g(s)\ < M\t — s\ a , —b<s, t < b. 


Let 

Then 


/ = Rg. 


\f(x)\ — » 0 as x — > b — . 


Proof 


We begin by noting that if 


then 


g i(x) = 1, |a;| < b , 

(Rgi)(x) = 


1 2 lb — x 


y/1 -M 2 r; \ b + x 


(vr) 


and — > 0 as x — »• b— . Now 


f b tti (#(£) ziMl 

l-b\b-Z Z-x 


df 


< M Jf fbTt 


\£-x 


\a—l 


Vb 


and 


/ je-^r 1 • (&-er l/2 de < 

for |rr| < 6, by a simple application of the Holder inequality, for 

<y — 1 = — T £, £ > 0. 

It follows that 


and hence so does 
as required. 


(Rg — Rgi)(x) — » 0— as x — > b— 
( Rg)(x ) — >• 0 as x — ► b— , 


(3.17) 
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Corollary 


( Rg)(x ) — > 0 as 

if g is absolutely continuous with derivative in L°°. Moreover under this condition, 

/ = Rg 

is the only solution of the integral equation: 

Hf = g 

such that /(•) is bounded at x = b— . 

It is convenient to rewrite (3.8) for a more detailed analysis of G(X,y). Since it is 
only a matter of scaling we will take 6=1, from now on. Let 

w(x,X) = w(x, A) 

i(;r, A) = A(x, A) exp 

K n (z) = zKi(z). 

Then (3.8) becomes: 

w(x,X) = j ^ A(Z, A) 

\ t , (X M \ 1 

| 2tt 111 (1—M 2 ) ^ ^7 (rr-0 

1 A / A M ..A 

+ -2nVT^W u Ao (t/ (1—M 2 ) 

~ WHF (^) X” exp (^rAe) *}' 

(3.18) 
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The crucial point of departure of our analysis is to rewrite the kernel in parenthesis 
in the integral above by 


Vl -M 2 1 


2n x — £ 


+ G'o(A; x — £) 


where, defining 


we have that 


^o(^; y) 


2n 


R 


11 


R n (z) = K n (z) - 1 
= zKi(z) — 1 


AM . A 1 1 

\y\ 


(3.19) 


A 


AM 


U (1 — M 2 ) |y 7 y ' 2kVT—M^ u K °{u(1-AP) 


1 


A 2 


27r7l— M 2 A 2 io 

Define the operator Go (A) by 


°° ( ~ Xt \ T ~ 

exp 1 tttz — rmr I A, 


AM 


U (1— M 2 ) / \U(1~M 2 ) 


|f - 2/| ) dt. 


(3.20) 


G'o(A)/ = g , ReA>0 

</(*) = G' 0 (A, rr— £) /(£) \x\<l. (3.21) 

4 _ 

Then Go (A) is readily verified to be a linear bounded operator on M into L°°, and 
(3.10) can be rewritten in operator form as: 

w(- , A) = H(A(-, A)) + G 0 (A)i(-,A). (3.22) 

What we have done is simply to isolate the “singular part” of the kernel as the first 
term 3). 

From (3.13) we can obtain our first breakdown of the solution by operating on 
both sides by R and obtaining 

Rw(- , A) = A(- , A) + RG 0 (A) A(- , A). (3.23) 
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Existence Theorem 


We can now prove existence of a solution of (3.23), which is actually an easy 

consequence of (3.23) where the main feature to be noted is that RGo(X) is compact 
4 - . 4 - 

on L 3 into L 3 . By the uniqueness of solution argument we see that 
(/ + RG 0 (X))A = 0 for some A in L 3 , 


would imply that 


HA + G 0 {X)A = 0 


or that (3.11) holds. But in that case we have seen that A must be zero. But 
RG 0 (X) being compact by the “Fredholm alternative” we have that (/ + RGo(X)) has 
a bounded inverse. Hence from (3.23) we have that 


A(- , A) = (/ + RG 0 (X))~ 1 Rw(- , X). (3.24) 


Or, we have both existence and uniqueness of solution. 


A Constructive Solution 

We can obtain a constructive solution for |A| small. Thus from (3.13) it can be 
seen that for each M < 1, 


||G' 0 (A)|| -> 0 as |A| -»• 0. (3.25) 

Hence for fixed M, 

||i?G'o(A)|| < 1 for |A| < |A m |. (3.26) 

Hence we have our first result on existence: 

Lemma 5 

For each M , 0 < M < 1, we can End |A M |, 0 < |A M | < oo, such that (3.15) has a 
unique solution given by 

A(- , A) = (/ + RG 0 (X)y 1 Rw(- , X) 

OO 

= Rw(- , A) + '^2(—l) k (RG 0 (X)) k Rw(- , A) (3.27) 

i 

for all | A | < |A m |. 
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Remark 1 

In particular we have from (3.27) that for A = 0: 


i(-,0) = Rw(-, 0). 


(3.28) 


Remark 2 

In [3] there is mention of an expansion given by Dietze [6] without proof. It is not 
clear to this author whether it is the same as (3.27) or not. Note that the expansion 
in (3.27) is not in powers of M . 

Remark 3 

Since the major interest in the use of the Possio solution is in the range of values 
of A such that 

k < 1 

where 


the expansion (3.27) is not totally void of value and we expect that taking the ex- 
pansion to the first term 

Ai (- , A) = Rw (- , A) - RGo(X) Rw (- , A) 

may be useful with further simplification of the kernel Gq(X) as well. 


Kutta-Joukowsky Condition 

We shall assume from now on that w ( ■ ,A), and hence w ( • ,A), has a bounded 
derivative in [—1, 1]. It would then follow from Lemma 4 that the function 

(Riu(- , X))x — » 0 as 

We shall now prove: 

Theorem 3.1 

4 _ 

Let / 6 L 3 . Then the function 

[RG 0 (X)f](x) — ► 0 as nl-. 



Proof 


Let us follow the breakdown of G'o(A) into the three terms as in (3.20). Let us 
also use the abbreviated notation: 

AM A 

7 “ U (1 — M 2 ) ’ " £/(l-M 2 ) ' 

Thus let 

(A) = G' 0; i(A) + G'o,2(A) + Go,3(A) 

with corresponding kernels 


G'o,i(A,y) — 

2tt 

. . 1 

li (7 2/ ) - 
y 

(3.29) 

G'o, 2 (A, y) 

1 

j; I<oh\y\) 

(3.30) 

2ttV1-M 2 

Go,3(A,y) = 

1 

\2 poo 

— 2 J 0 e-“ Kohlt -y\)dt. 

(3.31) 

27TV1-M 2 


Define now the operator 7 Zn(t) for each t > 0 on L "' by 


ft-ii (t)f = g(t , •)> f > 0 




v / 1 3 M I f 1 Ru(-y\x-Z-t\) 
2ir J- i x — f — t 


/(O M < i- 


(3,32) 


4 _ 

Then LLi(t) is linear bounded on LA into L°°. This is because the only singularity 
in the integrand is where 

x — f — t = 0 

but we can use the known expansion for K\ (z) at z = 0 and obtain: 
zKi(z) — 1 


A (2 Euler Gamma — 1) + | log - 


Re z > 0. 


(3,33) 


In particular we recognize that 


Go,i(A) = ftn(0) 

and that the corresponding kernel Go,i(A, y) is such that 

G 0 ,i(A, -y) = — G 0j i(A, y), y> 0 
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d 

Rnhy) 

d 

dy 

y 

dy 


lK \ {jy ) - - 

y 


y> o 


— + 7 

r 


iy 


- Kohy) 


y 


[1 - jyK^jy)] - 7 2 /v 0 (7I/) 


which at y = 0 using the expansion (3.33) 


= |(2 Euler Gamma — 1) + | log — — 7 s log 7 y 

which is in L p [— 1, 1] for p > 1 and hence it follows that 72-n(0)/ has an L 1 derivative 
in [—1, 1] and hence by Lemma 4 we have that the function 


[Rn n (o)f](x) 


0 as x — > 1 — . 


We shall abbreviate 7Ci(0) to Tin in what follows. We note also that for A > 0, using 
the expansion at y = 0 , 

Riiiiy) iRuiiy) 


7 2 y 


4(2 Euler Gamma — 1) + | log 


y 7 y 

we have the estimate, as a function of 7 that 

ll^n || = 0 [ 7 2 | log 7 | ], as 7 ^ 0 . 


(3.34) 


To handle the remaining two terms (3.30) and (3.31) we define the operator JC(t) 

4 _ 

on L * for each t > 0 by: 

K(t)f = g(tr) 


g(t,x) = 


1 


Kofrlt-x+tl) f(£) d£, |rr| < 1, 


(3.35) 


2-TrvG-Af 2 J- i 

4 _ 

Then JC(t) is a linear bounded operator on L 3 into L°°, and we have 

\ \ 2 poo 

G 0 (X)f = K n f + - j^J o e~ yt JC(t)f dt. 

Let us consider the second term. Let g = /C(0)/. Then integrating by parts, we 
have: 


g( x ) = 


1 


2tt Cl -M 2 


1 K 0 ( 7 (l-x))F(l) - 1 1 


m 




+ 


(1 -M 2 ) 


2n Cl -M 2 J-i x - £ 

(Tin F) (x), |x| < 1, (3.36) 
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where 


We define a new function A" (7, •), for each 7, by 


K{n,x) 


1 /1 /i /rrr A: 0 ( 7 (i-e)) 


7r 2 V 1 + x 7-i y 1 — £ £ — a; 


|x| < 1 . 


Then it is immediate that 


R 2 tt /T^M 2 Ao(7(1 (1 —M 2 ) A ^ 7, 


(3.37) 


Lemma 5 


K (7, a?) — s- 1 as x * 1 


(3 .38) 


Proof 

tfo(7(l - 0) + log(7(l -0) 
is Liptschitz in |£| < 1. Hence 

/1-s r 1 FT (JfoWi - 0) + iog7(i - 0) 


1 + x 7-i y 1 + £ 




and so does 


/ 1 ~g f 1 1 1 7 

\ 1 • .r 7 . v 1 + e £ .r 7 

By the Corollary to Lemma 3 we see that 

K (7,21) — > 1 as x — + 1— . 


0 as 


Hence 


7? [A(0) /] 


_ £(o , 1 o r7? 

(1 —M 2 ) (1—M 2 ) (1— M 2 ) 11 J ' 


(3.39) 


where every term is bounded at 1 — , as required. 

Finally we consider the third term in (3.31). The corresponding operator is given 

by 

— \2 /.co 

Go, 3(A)/ = — Y J o e- xt JC(t)f dt. 
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Let 


g(t,-) = 


Then integration by parts yields: 


, ^ 1 1 
9{ti x ) = o 7= 


2tt Jl-M 2 


K 0 ^{l+t-x)) F(l) - — 


1 1 /•' F(f) 


27T y/l—AT 2 J-i x — t — £ 


1 1 
7T \/l-d/ 2 


(Fn(t)F) (x), |x| < 1, 


(3.40) 


where 


F(x) = 


We carry out the integration with respect to t, term by term. Let gi(t,x) denote the 
first term in (3.40). Then 


-F f e->‘ 9l (t, X ) dt = -±± fe-*'KoMl + t- X )) dt 


-A 2 1 F(l) r°° _ } 


which, by integration by parts, 




e xt r fKi('y(l+t—x)) dt 


A F(l) 

U 2tt 


T -Ad 2 K 0 ( 7(1 -x)) 


U 2tt ] do 


130 w, dt 

e~ A * 


1 + t — x 


- I< 0(7(1 -x)) 


+ / e' 


-Vt ^11(7(1+^— 
1 + t — x 


dt) . (3.41) 


The third term in (3.41) 


77 V 1 — M 2 
U 2tt do 


F(l) d°° —x't RiMl+t-x)) 


1 + t — x 


dt, |x| < 1 




Rail 1 ) 


4 V / 1 Z M 2 e A'(i-a:) r e -t/M (U 

U 2?T J 7(1— x) 


where 


Fi(t) 


|t|Fi(|t|) - 1 
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Hence defining 


'i-M 2 r 00 


,X(1— x) 


roo 

/ e -*/ M m(t) dt , |x| < l 

J 7(1—0;) 


the third term can finally be expressed as: 




/*CO /‘OO rj 

/ e~ t/M r n (t) dt = / e~ t/M — [-I< 0 (t) - log t] dt 

Jo Jo at 

1 /»oo 

= -yy y o e“ t/A/ [-#„(*) - log/:] dt 

roo 

= / e -t [— Ao(Mt) - log M£] dt 

Jo 


(- 1 ) , 

>\ — M 2 Gg M M 


log M + log 2 


=i= (l-Vl-M 2 )logy - log 


and is < 0 for 0 < M <1. Hence 


roo roo 

/ e~^ M rn(t) dt = / e _ fi M (— ?’n(f)) dt 

Jo Jo 


1 T (l + V 7 ! 3 ^) ( r - — — — \ M 

7 log — ( 1 — \J\ — M 2 ) log — 

r \-M 2 & 2 v ) & 2 


_ "M 2 1 1 

O log as M — *■ 0 

2 & 2 M 


as M 1. 
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Hence for 7 > 0 (equivalently A > 0 ) 

Af( 1 ) 9 3 W < Jj \F(l)\e 2x ' e~ ,,M (~r u (t)) dt 

A M 2 1 

~ j 7 ^|e(l)| — log— as M —> 0 (3.43) 

~ F exp ( w(klFj ) |F(1)I “ "- 1 - (3 ' 44) 

Note also that the function q 3 (') does not depend on /(•)• It is also immediate that 

Rq 3 ) (x) = 0, xml—. 



where 


Hence 



(3.45) 
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is given by the function: 


i r p m 


K( j,x) + - \ J—~ h( x ',x) 

7T V 1 + X 


+ A F(l)(Rq 3 (-))(x) 


Lemma 6 
Let 

Then 


h( A, x) 


1 1 — x 


IT V 1 + x 


h(X',x), |x| < 1. 


h{ A, x) 


— 1 as x — > 1 — . 


Proof 


1 1 — x 


h(X,x) = / 

7T y 1 — X Jo 


TT \ 1 + X 


-At 


/2 + 1 


t + 1 — x 


dt 


— 1 . ■ re-^-^J2 + P(l-x) 

7 r VI + z Jo v 


VT + 

— 1 as x — >■ 1 — . 


Combining with Lemma 5 we obtain: 

Corollary 

The function defined by (3.46) goes to zero as x — > 1 — . 

Next let us tackle the second term in (3.40). Here we note first that 


F(0 




1 -t Fin r 1 

v ' df + 1 


m 


J-i x — t — f J- 1 x — t — f 

The first term on the right can be expressed: 

r 1 nj - 1 ) 

Jt-l (2: ~ C.) 


1 -t x — t — f 


df- 


df 


(3.46) 


(3.47) 


dt 

1 +t 2 


can write: 
(3.48) 
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which in turn can be expressed as 


q(t,x) 


F(Z - 1) 


defining 

Hence 


i (x-£) 

m = o, iei > i. 




R 


1 1 


?(V) 


(x) 


F(x — t) 
(1 — M 2 ) 


_2t r V / T Z M 2 

Let r(t,x) denote the second term in (3.48). Then we can write 

d F(l — £+a) 


•(M‘) = / 

Jo 


x — 1 — (J 


da 


and 

1 1 


2tt vH-M2 




1 1 2 1-x 


2n \/\—M 2 7T V 1+z do 



A. / 1+ « 1 


i V i~ £ 1— <j 


1 1 l-z d 


1 


/ 2 + <7 


1 — M 2 7T V 1+# do x—l — cr V CT 


F(1 — £ + <t) da. 


Hence 


-A 2 1 


,-A't 


H 2 

is given by the function 


i? 


- r 

2 7T 7-1 


f«) 






A 2 


1 


U 2 1 - Ad 2 do 
A 2 1 


e * F(x — t) dt 


. . e V *F(1 - t) dt] h(X',x) — \i 

f/2 1 - M 2 Vdo V ' d V ’ 7T V l+x 


which upon integration by parts in both integrals 

-F(x) + r e*C~> f(a) da 


-A 

IT 


- jj h(X,x) 


~F( 1 ) 


^Al- 


and clearly goes to zero as x — > 1— . 


F(1 — £ + a) d£ 


dt 


x\ < 1 (3.49) 


/(o') da 
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Finally let us consider the third term in (3.40) and calculate 

-A 2 


R 


roo 

/ e~ xt R 11 (t)F dt 
Jo 


_U 2 (1—M 2 ) 

Let g.(-) denote the function inside square brackets: 

1 1 no a 2 r u t Rnh\t-(x-o\) 


n x ) = / br 


1 2 tt \J\—M 2 
1 1 A 2 F(0 




where, 


= /_>-( 


i 2 tt U 2 ^l^M 2 

llARiillA) “ 1 


t~ (x ~0 
L(x - 0 e~ x '^ x ~^ d£ 


dt , |x| < 1 


,-A 't 


0\t\i<i(0t\) - i' 


dt 


dt -f- 


mi 


m KiOf\t\) - r 


dt 


where the first integral, denoted c(M ), is given by (3.42). We can decompose g 4 (-) 
correspondingly 

? 4 ( ;:r ) = J-i.iO) + (3.50) 


where 


? 4 ,lW 


7 4 ,2^J 


1 c(M) v , ( A 2 /■ 

2 tt VI -M 2 


m f PM'*# 


(3.51) 


1 A 2 


-At 


1 


2? r t/ 2 Vl-M 2 

• J' e*F(OdS j 


-X't 


'7l*|Ki(7l*D - 0 


0-£) 


7 


log 7 


as 7 — > 0 . 


dt (3.52) 
(3.53) 


It is immediate that 


{Rq 41 )(x) — t 0 as nl- 

(. Rq 42 )(x ) — >• 0 as x — > 1— . 


Hence, finally: 


i? 


-A 2 

H 2 


e A * KL(t)f dt 
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is given by the function 

-^F(l)[l<( 7,ar) + h( A',x)j - ^ -F(x) + J^e x ' (x - a) f(a)da 
- jj ~F( 1) + J^ e - x '( x -<r) f(a) da h( X',x) + (Rq 4 (-))(x) 

= -jj- .F(l) K ( 7 , x) - j - F(x ) + j ^e _A ' (x_<r) f {a) da 
~ jj f x e _A(1_<T) f(a) da h{ \',x) + (Rq 4 (-))(x). 

Hence 

A A 2 r°° 

RG 0 (X)f = RKuf + R jjF(0)f - —J o e-X dt (3.54) 

where 

\ \2 7*00 . 

R V mf - JH l e- xt IC(t)fdt 

is given by the function 

[K(l,x)F(l) - F(x)] - jF(l)K(^,x) 

-j[~F(x) + J\- X '^f(a)da 

- Tr (/ 1 e_A(1_ff) /( 0 ') da ) ~h(x',x) + (Rq u (‘))(x) + F{l)(Rr n (-))(x) 

= ^ /V*c~> m da 

-jh(X',x) J ^ e“ A ' (1_(r) f(cr) da 
X M 2 

+ uo^m [F{1)KM ~ F{x)] 

+ (Rq 4 (-))( x ) + JJ F(l)(Rq 3 (-))(x), |x| < 1. (3.55) 

Using Lemma 4 (or otherwise) we see that the function defined by (3.55) goes actually 
to zero as x — *■ 1— . This concludes the proof of the Theorem 3.1. □ 
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Closed Form Solution Approximation 

We now show how we can derive a closed form solution with error of the order 

|') 2 log 7 | as 7 — >• 0 


\M 2 log M\ as M^O. 

For this purpose we note that by (3.55) we can break RGo(X) as 

RG 0 (\)f = T(X)f + (RKu(0)f + Rq 4)2 ) (3.56) 

where the part in parentheses, as we have seen, 

0 [ I 7 2 log 7 I ] as 7 -H- 0 

and T( A) is defined by 

9 = T(X)f 

g(x) = I ^ e - A ' (;r ” (T) f(a) da - ~ h(X',x) j ^ e “ V(1 “ CT) f(a) da 

X M 2 A 

+ jj l _ Ar2 [i r (l)/A(7,x) - F(x)} + — F(l)(Rq 3 )(x) 


+ (Rq 4il ) (x) 

where the last term can be simplified further, noting that 


jG' s F(t)d( = A (e v F(l) - /V £ /«)*;) 


so that 


(Rq 4 , i)(x) = c(M) ~ F{l)q 5 (x) - q 5 (x ) j^e A(1 0 /(0 


where 


9s( 


1 1 — x 


x) 


l+c e 




7T 2 v 1 + x \.l 1 y 1 — £ £ — x 

both of which are special functions, and do not depend on /(■). 




(3.57) 


(3.58) 


(3.59) 
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Collecting common terms, we can rewrite T(X)f as 


g(x) = f t e X ' {x a) f(a) da - 7 M J j{a) da 

— 7 T h(x) f e - A '( 1 -°') J((j) da + r yMK(x)F(l) : |x| < 1 (3.60) 

U 7-i 

where 

h(x) = h( X',x) + c(M) q 5 (x) 

Y 47 2 

K{x) = K( j,x) + M2 ((Rq 3 )(x) + c(M)q 5 {: r)) . 

H(-) and K(-) are in L 3 . 

The main point in isolating T( X) is to exploit the fact that we can obtain a closed 
form solution for 

(1 + T( A))- 1 . 

Since the additional terms in (3.56) are small compared to T( X) for small 7 , (or small 
M for fixed A), we obtain in this manner a closed form solution for the Possio equation 
which is particularly suited to the usual range of 7 values of interest: I 7 I < 1 . 

The reason why we can get a closed form solution for g for given / for the equation 


/ = (. I + T(\))g 


is that (/ + T(X)) has the form: 

(I + T(X))f = g + Lg + ftC^g) + f 2 C 2 (g) (3.61) 

where 

Lf = -jj- J ^e~ x ' {x ~ a) f {a) da - 7 M J J (a) da (3.62) 

is a Volterra operator, £i( - ), £ 2 ( - ) are gi yen linear functionals and 

fi = h, f 2 = K 

4 

are given elements in L 3 , and we can evaluate (/ + L) 1 in closed form. Hence 
9 = (I + L) + (I + L) + (I + L) ^^(g) 

and we only need to solve an algebraic linear equation in two unknowns for C\ (g) and 
£2 {9)1 which we can do provided of course the determinant is nonzero. 

It is a little easier to proceed directly as we shall show. 
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Theorem 3.2 


The equation 

f = g + 7 (A)// 

has the unique solution given by 

/ x rx 

cosh 7 (x — a)r(a) da + 7 J sinli 7 (x — a)r(a) da (3.63) 

where r(x) is given by 


r(x) = f(x) + 7 S(r)K(x ) + (— h(x) — 7 MK(x)j C(r) 


where the coefficients (functionals) 

S (r) = j sinh 7 (l — a)r(a) da 
C(r) = J cosh 7(1 — a)r(a) da 

are given by 

I C(r) I 

= d(\ My 1 


(3.64) 


(3.65) 


C(r) 
S(r ) 


<?(/) 

5 (/) 


(3.66) 


where the 2 x 2 matrix 


D{X,M) 


1 - A C'(h) + 7 MC{K ) 7 C(K) 


~S(K) 


1 + 'yS(K) 


(3.67) 


is nonsingular for Re A > 0 excepting possibly a sequence of isolated values { A^ } 
bounded away from zero and | A* | — » 00 as k — > 00 . 

Proof 

We begin with 

f(x) m g(x) - J ^ e" V( - T ~ <7) g(a) da - ^ h(x) J ^e~ x ' (1 ~ a) g(a) da 


+ M 2 A' 


rl rx 

K(x) j g(a) da — J g(a) da 


(3.68) 
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Let 


/>) = e x ' x f(x) 
g(x) = e x ' x g(x) 

/ X 

9{ a ) da. 

Multiplying both sides of (3.68), by e Xx we have 

f(x) = g(x) - ^ j J(a) da - e x ^ x ~ 1] h(x) J^g(a) da 


+ MX 


Now 


K(x) J e x '^ x ^ g (a) da — J e x '^ x ^ g (a) da 

J[ ‘ e xl( - x - a) g (a) da = Q(x) + X' j X e V(a:_<7) g(x) da 
J 1 e x '( x -*) g ( a ) da = Q(l)e A 'M + X f ^ x ~ ff) Q(a) da. 


Hence we can write 

f(x) = Q\x) - ^ Q(x) - Q{1) ^ c V(a '- 1) h(X\x) 


+ M 2 X e x '^ x - l) Q{l)I\{x) + X 2 M 2 f e x '( x ~^Q(a) da I<( 
- XM 2 Q(x ) - X 2 M 2 J X e x,{ - x - a) g(a) da 

Q'(x) - X'Q(x) + Q( l)e A ' (x_1) 

K{x) ^ e x '^- ff) Q{a) da - J* e x ' (x ~ a) Q(a) da 


XM 2 K(x ) - A h{x) 


+ X' 2 M 2 


Q'(x) = XQ(x) + e x ' x f(x) 

+ g(i) e M^- 1 ) 

+ M 2 X' 2 


p h(x ) — \'M 2 K(x) 




Q(a) da — K(x) J e A ( - x Q(a) da 
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which yields 


Q(x) = J ^e x ' x f(a) da 

+ M 2 A'" J e X x da j e~ x ' T Q{t) dr — K(a)e x ’ x j e~ x ' T Q(r) dr 


e x ' x e V — h(a) — \'M 2 K(a) da. 


Hence letting 


we have: 


Q(x) = e x ' x Q{x) 


Q(x) = J J(a) da 


/*CC c 1 ™ 

(x — r)Q(r)dr — J K(a) da j Q(r) dr 


+ Q(i) 


(— h(a) — X 'M 2 K(a) \ da 


where we have used: 


J da J e X ' T Q(r) dr = J (x — r) Q(r) dr. 


Hence 


Q(x) — X ' M 2 J (x — t)Q(t) dr 


J f(a) da — X' 2 M 2 (^ j K{a) da^j (^ j Q{r) dr 
+ <5(1) / ( 77 H a ) ~ X' M 2 K (a)\ da. 


This is a Volterra integral equation which has a unique solution. In fact let 


Then 


F{- 1) = 0 = F'(-l) 
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and 


F"(x) - X' 2 M 2 F(x) = r(x) 

r(x) = f(x) - X , 2 M 2 K(x)F(l ) + F'(l)(\h(X',x) - X'M 2 K{x)) 


yields 


F(x) 


x sinh 7(2; — a) 


7 


r(a) da , 7 = MX 1 


or, 


F'(x) = I cosh ^{x — a)r(a) da 


rx 

Q(x) = J cosh 7 (x — a)r(a) da 

I f X 

Q(x) = COSI17 (x — a)r(a) da. 


Hence it follows that 


/ X 

cosh 7(7: — a)r(a) da 

/ X 

sinh 'f{x — a)r(a) da 


where 


r(x = 


f(x ) — 7 K(x) J sinh 7(1 — a)r(a) da 
+ h(x) — 7 MK(x)^j J cosh 7(1 — a)r(a) da 


where we need to determine the coefficients C(r), S(r), where C(-), S(-) 

4 _ 

tionals defined on L*- by: 

C(f) = 


ar 


S(f) 


cosh 7(1 — a)f(a) da 


sinh 7(1 — a)f(a) da. 


But from (3.70) we obtain the algebraic linear equations 

D(X, M) 


C(r) 


C(f) 

S(r) 


S(f) 


(3.69) 


(3.70) 

func- 
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where D( A, M ) is as defined by (3.67), and has a unique solution provided 


d(X,M) = det D(A, M) 

= (l - ^ C(h) + 7 MC(I < (1 + 7 S(I<)) + 7 C(K)S{K ) 

= (l - ^ C-w) + 7W (i - ^ C-wl (3.71) 

+ S(K)C(K)( 7 + 7 2 M) 


is nonzero. 

We note the d(\,M) is continuous in Re A > 0, and analytic (no poles or other 
singularities in the finite part of the plane) in Re A > 0. Moreover 


d(M, A) — > oo as ReA ^ oo 


and hence it has a finite number of zeros in any finite part of Re A > 0, the sequence 
{ Afc } of zeros being such that | A* | — > oo as k — > oo. It is a continuous function of 7 
in (3.71) and for 7 = 6 (or M = 0): 


= 1 




H A', •) 


= 1 - uL h {i M da 

1 

Ci (A) 


POO 

Ci (A) = / e~ xt Ci(t) dt 
Jo 


U e- x ' u 

1 K o (if) + K, (A) 


and 


ci(A) — > 00 as Re A— > 00 . 

Hence it follows that the { A* } are bounded away from zero. This concludes the proof. 

□ 


Remark 1 

We conjecture that d( A, M ) has no zeros in any finite part of Re A > 0. 
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Remark 2 

Let us summarize our main result. We have derived an approximate closed- form 
solution to (3.8). 


where we set 


A(x,X) = exp( ^^^ jA(x,X) 
A(x, A) = g{x) in (3.69) 


f( x ) 


w(x, A) exp 


/ -A M 2 x \ 
\U(1-M 2 ) ) ' 


Remark 3 

For M = 0 (incompressible case) we have that 

T( A) = RGo(X) 

and further (3.42), (3.43) reduce to 

A f x 

g(x) = r(x) + — via) da 
U J - 1 

r(x) = f(x) + ^ h(j^ , x^J J^r{a) da 

and 

J r(a) da = c x (A) j f(a) da. 

Hence 

g(x) = r(x) + — r(a) da (3.72) 

U J - i 

where 

r(x) = f(x) + ^ c x (A) J J (a) da. 

This checks with the known result for the incompressible case, as we shall see in the 
next section. 

Remark 4 

To simplify calculation without losing too much accuracy we may take 


= K{ l,x) 

= h(X',x). 


K (x) 
h(x) 
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Velocity Potential 

The Laplace transform of the velocity potential, (j)(x,z,\), satisfies 

8 ~ 

\(/>(x,z,\) — (j>(x,z, 0) + U—(f>(x,z, A) = i/;(x,z, A). (3.73) 

ox 

Solving this differential equation for <j)(x,z. A), we obtain 


A) = exp 


— X(x — a ) 
U 


4>{a,z, A) + 


l r x 


U Ja 


exp 


— A(a; — s ) 
U 


^o(s, -2, A) ds 


where a is arbitrary. From (3.3) we see that 

lim <p(a, z, A) = 0. 

a— oo 

Hence 


4>(x,z,\) h — J 
Hence in particular for z = 0+, 


exp 


— X(x — s ) 


U 


( x , 0 


A) = -r£ exp ( 


— A(.t — 5) 
U 


ip 0 (s, 2 , A) ds. 


A{s, A) ds. 


(3.74) 


Also, differentiating with respect to x: 


- 1 ) 


d<p(x, 0, A) 1 / — A N 


dx 


2 V u 


x ( —\{x — sY 

- b exp \— r - 


A(s,\)ds + b A(x, A). (3.75) 
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4. Subsonic Incompressible Case 


We now specialize the Linearized Compressible Flow equations (2.1) to the case 
where all the time derivatives therein are set equal to zero, yielding (by cancelling 
out cioo , it being nonzero) 


(1 — M 2 ) 


d 2 4> d 2 (f) 


dx 2 


+ 


dz 2 


(4.1) 


but (/>{•) is still a function of time, since the boundary conditions are retained: viz., 
the flow tangency condition, which then makes the boundary value a function of time 
through the downwash condition, and the Kutta-Joukowski condition, which also 
depends on the time coordinate, since the acceleration potential is still defined by 


ip(t,x,z) = 


d<t>(t,x,z) d<j>(t,x,z) 

dt dx 


However (4.1) is not an initial value problem. We are satisfied with any solution. Note 
that the velocity U has no meaning to (4.1), and enters only in the Kutta-Joukowski 
condition. 

This problem with M = 0 is the celebrated “incompressible inviscid flow problem.” 
(See [3,8].) We need to consider it for M / 0, and this is traditionally done through 
the so-called Prandtl-Glauert coordinate transformation on the truly “incompressible” 
solution for M = 0. (4.1) is also referred to as the “steady state” equation since time 
is not present. However it has nothing to do with “steady state” in any sense of 
asymptotic behavior in time. 

The main feature for us is that we can indeed obtain a “time domain” solution 
unlike the Linearized Compressible Flow equations where we can only obtain the 
Laplace transform of the solution. Even though no claim is made that the derivation is 
new, we shall need to be precise mathematically (as compared, say, to the “standard” 
treatment in [3]). We shall also present a few results which do not appear in [3] or 
elsewhere. 

We shall present a direct approach without recourse to the Prandtl-Glauert trans- 
formation. Thus, we shall simply begin, as usual, with the general form of the solution 
to the “potential equation” (4.1): 


(f>(t, x, z) 


-1 1 
27t v 7 1 —M 2 



7a (b 0 tan 


zVl-M 2 \ 

( x ~0 ) 




1 , fb+Ut ( Z J 1 —M 2 \ 

- 2 ? ' /TQ ^ L 7 » <( ’ e) tan " (iber) ^ <42) 

where the functions 7 „(£,■) and 7 are to be determined from the boundary 
conditions. 

From (4.2) we readily deduce that 
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d(f)(t,x, 0+) 
dx 


0, 

0, 


x < — b 
x>b + Ut 


\% a(t,x), \x\<b 


= |t w (t,x), b<x<b + Ut. 


Define (the “circulation” ) : 


r (t) = [ To(D x) dx. 
J-b 


Then, superdots indicating time derivative, 

d r b dcj)(t,x, o+) 




dx = (p(t,b—, 0). 


l-b dx 

As usual, we now determined Tu>(L •) from Kutta-Joukowski conditions. Let 

fb+c d</>(t,X, 0 + ) 


F(t,0 = -2 1 

Then we recognize that 


dx 


dx , — oo < £ < oo. 


fb+C 

F(t,() = J T w(t-,x) dx for ( > 0. 


The Kutta-Joukowski condition 

dcp(t , x , 0) 


yields 


0 


d 

dt 


dt 

b +c d(j)(t,x, 0) 
b dx 




dx 




+ Uh 7u,(t, 6+C), C > 0 


or 


or, 


|f(t) + - — F(t, C) + T'Ty(Lfr+C)> C > 0? 


0 = jfft) + <>o. 


or, we have the partial differential equation 

9F(t, C) _ v dF(t,0 


dt 


‘X 


- r(t), !>«,(>«. 


(4,3) 


(4.4) 


(4,5) 


(4,6) 


(4,7) 
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We note that 


F(t,() = F(0,(-Ut) + / T(a)da (4.8) 

Jo 

is a solution of his equation for ( > 0, and by uniqueness of solution of the Cauchy 
problem associated with this equation, the only solutions such that 


In particular therefore 


and hence for ( > 0: 


[b+C 

F( 0,C) = / 7u>(0 ,x) dx, C > 0. 

Jb 

dF(t, Q = dF( 0, C~ Ut) 
d( d( 


dF( 0, C- Ut) 


-- 1 w{t, b + Q. 


Let us use the notation 


5(0 : 


9F{ 0,0 


-oo < C, < oo. 


Then 

7 w(t,b + () = g(C~Ut) 

and going back now to the Kutta-Joukowski condition (4.6), for C = 0, therein: 

t(t) = -Ug(-Ul), I > 0. 

But 

7 w(t, b + C) = g( C - Ut), C > 0 


7 w (t,x) = g(x — b — Ut), x > b. 

Hence for 

x — b — Ut < 0 or x < b + Ut, 
we have from (4.9) that 


f u TT+\ f> fUt + b + X 

g(x — b — Ut) = — T I 


— 1 • / x — b\ x — b 

TT r< -— ’ ■ 


1 fw(t,x) = 


=±r(t - X x± 
v \ u 


b < x < b + Ut 


(4.10) 
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which is then one of the primary results of this theory. 

We note that the Kutta Condition at b— requires that 

f(t) + U la {t,b -) cfg 0 (4.10a) 


since by (4.3): 


d(j)(t, b— , 0 +) 8 f b ~ d(/)(t.,x, 0 +) 


dx 


r (t) 


dt dt J-b dx 

and (4.10a) is consistent with (4.10), upon taking x = 6 +. In particular we have 

n fw(t,b+) = 7 a (t,b-). 

Next we exploit the flow tangency condition. We have 
d(j)(t,x } 0 +) 


w(t,x ) = 


dz 


* (zi f d( _ ' f >+m idMl ^ 

\2k J-b X — E, Z7T Jb X — t; y 

where the second term, substituting for 7 w (t-,£) from (4.10): 

Z7T ,/6 U X — t, 


Hence 


re(f, a:) 


1 [ b 7a(t,0 


-f 

2n J- 


d£ + tt: 


2ty Jo x — b — Ua 
1 d T(t — a) 


da 


(4.11) 


\/l — M 2 2ir J-b x—l; 2tt Jo x — b — Ua 

which is then the “airfoil” equation, a singular integral equation which we need to 
solve to determine 7 a (t, £). We have covered the necessary relevant theory in Section 3. 
Rewriting (4.11) as 




7 a(U) 

-b X-£ 




tJ , > 


—M 2 -l— da — w(t,x). (4.12) 

x — b — Ua 


We have, under the condition of Lipschitzianness, and introducing the operator R , 
defined by (3.14), for the only solution bounded at £ = b— : 


7 a(V) = -Rw(t,-) + q(t, •) 


(4.13) 
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where we define 


q(t,x) = 


7r 2 V b + x J-b v b — £ £ — x Jo £ — b — Ua 


1 

jb-x 

7T 2 

V b + x 

1 

jb — x 


1 


ft T(t- a) 


da d£ 


2 + U U , 

— w- / — ai— dcr, U = — . (4.14) 

7 t b x Jo x — b — Ua ^ Ua b 

Remark 

In the aeroelastic problem (as in [3]), the downwash function 

w(t,x) = — h(t) — (bx — a) a(t) — Ua(t ) 

where /?.(•) is the “plunge” and «(•) the “pitch” angle. 

From (4.13) setting f = 0 we obtain 

7o(0,-) = ~Rw( 0, •) (4.15) 

showing that the initial aerodynamic flow conditions are determined by the initial 
conditions of the downwash. Our next step is to determine the circulation function 
from (4.13). It is convenient to let 


and 


so that 


B(t , •) = — Rw(t ,•) 

„ r b 

B(t ) = J B(t,x) dx , 


r(0) = B(0) 


and, integrating (4.13), we obtain 


r (t) = Bit ) + / q(t,x) dx 

-b 


where 


r b ft 

<l(t,x) dx - / 

-b Jo 


1 - 


2 + Ua 

N Ua 


or, 


0 = 5(f) - B{ 0) - f 

Jo 


2 + Ua 
N Ua 


r (t — a) da 


T{t — a) da 
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which is a linear integral equation for T (t), where the important observation is that 


A / e 
J o 


2 \-Ut 


-At ^ ^ ^ 

\ Ut 


dt 


ci (A) 


where 

ci (A) 
and hence 

Af(A) - r(0) 


u 

e-x/u 

= r 

x K 0 \ 

4) + K l (4) 
\uJ 1 \uJ 

Jo 


roo 

/ e~ xt Ci(t) dt , Re A > 0 

Jo 


(4.16) 


B( A) 


1 


B(0) ' 

A ) h( A) 


: ci (A) (\B(X) - 5(0)) (4.17) 


yielding a “convolution” interpretation in the time domain. This is another central 
result in the airfoil theory, the key being (4.16) due to Sears (see [3]). However this is 
not of particular importance in the Laplace transform theory for (4.13). Thus Laplace 
transforming (4.13) we have, using (4.17) as an intermediate step: 


7 a (A, •) = -Rw( A,-) + (AT(A)-T(O)) h(\-) 


or, using (4.17), 


7 a (A, ') = -Rw( A,.) + (AB(A)-B(O)) c,(A)ft(A,.) 


(4.18) 


where 


h( A, x) 


1 jb — x f°° 

1 

2 + Ua 

iryb + x Jo 

x — b — Ua \ 

POO 

Ua 

7a(A,x) = j 

ss 

TT 

-f-T 

cS 

■w 

< 

1 

03 

O 



POO 


B(A) = j 

1 e~ M B(t ) dt 
0 



da 


To finish the calculations we need next to evaluate 7 w (t, •), using (4.10). We have 

7 »M + C) « ^ r^- o<c <ut. 

But the time domain version of (4.17) yields 

T (t) = [ Ci(t — a) B (a) da 
Jo 
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and hence 


-1 ft-«/u) 

7 w (t, b + Q = — ci(t — a) B(a) da 

U Jo 


also 


o, c > ut, 


7 w (A, b + Q = e + w (t, b + Q dt, 0 < ( < Ut 
Jo 


-jr e-* ,U (Af (A) - r( 0 )) 


u 

= ^e- A " c 'c,(A)(AB(A) - B(0)), 

This yields for the acceleration potential 

■0(A, x, 0+) = %{X,x) + J b (X^ a (X,y) - 7a(0,y)) dy (4.19) 

using 

%(X,x) = (-Rw( A, •))(*) + (XB(X) - B( 0)) ci(A) h(X,x), (4.20) 

and (4.15) 

7a(0,y) = -(Rw(0, -))(y). 

We note that (3.49) is the same as (4.19) if we set the initial conditions 

7a(0 ,y) = 0; B( 0) = 0 

and take 

r(x) = 7 „(A,.c) 

/(•) = -Rw( A,-) 

and in fact this analysis provides an interpretation for the function r(-) in (3.72). 
Also, finally, we have for the velocity potential: 


(f)(t, x, z) 


-1 f b , ^ z\J 1 —M 2 „ 

~r~ / 5(i,^)tan — — 

Z7T 7-6 [X — Q 


1 /* 1 /•* r(t-u) 


2 + f/u _1 zQT^KP 

‘2tv J-b n\l b + £ Jo £ — b — Ua \ [Ja (x — £) 


1 A • x Wl-Tf 2 

+ — / r (£ - cr) tan 1 — — da, 

27 r 7o (x — 


(4.21) 
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where 


T(t) = [ C\[t — a) B(a) da 
Jo 

B(t , •) = —Rw(t,-) 


B(t) 

e~ xt ci(t.) dt 


B(t,x) dx 


ci (A) = 

J U 

Note that the initial conditions are: 


oo TJ 

At - ' - u 


,-A /U 


X K ° (f) + K ' (f) ' 


_ 1 r b , 2 a /1 — M 2 

<^(0 ,x,z) = ^(0,0 tan- — - df. 


0(0, of, 2 :) = 2^ J_ b mo tan 


Hence 

On the other hand 


(®-0 

.! zVT-m* 
0~0 

<t>( 0, rc, 0) = 0 = c/)(0,x,0). 


dO 


^(0+,ar,0+) = 7„(0+,:r) = B( 0,x) 

and is not necessarily zero. Hence 

' dd> dd> 

°’ x ’ 0) - £, a; + m 

= UB{ 0,x) = 0, if w a (0,-) = 0. 

To connect with the Possio solution (cf. (3.49)), we note that 

ip p (0,x,0) = i/j o (0,x,0) + cf) p ( 0,x,0) 

= A(x, 0) = 0, if w a (0, •) = 0 = w(0, •)• 


Hence ^/’(O, rc, 0) and i[’ p (0,x,0) are equal if the downwash initial conditions are zero, 
but not otherwise. 
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5. Solving the Boundary- Value Initial-Value Problem 


Solving the Boundary- Value Problem 


We first consider the problem of finding a particular solution to the Boundary 
Value Problem for (2.1) without regard to the initial value. Thus let A L denote the 
extension of A where we drop (only) the condition (i) that the boundary value be 
zero. Thus 




fi 


9i 

AlS = 9 

f = 

/a 

9 = 

92 



h 


93 


where 


-2 U& al(l-M^ x 

0 0 

0 0 


JL 

dx 


d_ 

dz 


/■ 


(5.1) 


The problem then is: Given w(t,-) where for each t > 0, w(t, •) £ L 2 DR 1 ], find f(t ) 
in 7i such that 

i) A L f(t) = 0 


and 


ii 


/ OO 

\fs(t,x,z) — w a (t,x)\ 2 dx — >■ 0, as z 0+ 

-OO 

It is easy to verify that such a solution is given by 


(5.2) 


f(t) 


/i(-) 

/a(-) 

fs(') 


fi(t,x,z) = 0 


d 


f 2 {t,X,z) = —<f>(t,x,z,) 
d 

h(t,x,z) = —<j>(t,x,z,) 
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where <^(t, x, z ) is the velocity potential determined explicitly in Section 4, equation 
(4.21), satisfying (4.1) and the flow tangency condition and the far field conditions. 
We shall denote this velocity potential by < ) where M is the subscript to 
indicate that M is not taken to be zero. We can also calculate the velocity potential 
via the Possio equation as in Section 3, setting M = 0 in the acceleration potential 
(but retaining M in the definition of R). As we have seen, the two determinations 
differ in the initial conditions. In what follows we may pick either one and still denote 
it 4>m{')- 

We note while f (t) is not in the domain of A$ , it does satisfy the dynamic Kutta- 
Joukowski condition by our construction. In other words 

f x d 

'*l’{t,x,z) = J — f 2 (t, y, z) dy + Uf 2 (t,x,z ) 

will satisfy 

1p(t; 6-,0+) = 0 
ip(t, x, 0+) =0, x > b. 

Moreover we can write 

f(t) = Dw(t, •) 

where D defines a linear bounded transformation on L 2 \R i \ into 7i for each t > 0. 

Solving the Combined Problem 

We can construct the solution for the problem (2.1) — the combined initial- 
value/boundary- value problem, following the general technique outlined in [1,2,3]. 
We claim that the solution is given by defining: 

x(t) = x 0 (t) + Dw(t,-), t> 0 (5.3) 

where 

x 0 (t) = S(t)(x( 0) — Dw( 0, •)) — [ S(t — a)Dw (a, ■) da (5-4) 

Jo 

(the superdot as usual denoting time derivative) which is the generalized (or weak) 
solution of 

Xq ( t) = Ax 0 (t) — Dw(t , •) (5.5) 

(with the tacit asumption on strong differentiability of Dw(t , •)) with the initial con- 
dition 

Xo(0) = x(0) — Dw( 0,-). 
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To verify that x(t) given by (5.3) is the solution sought, we have only to note that 


x(t) = ± 0 (t) + Dw(t,-) 

= Ax 0 (t) 

= A L x 0 (t ) 

= A L [x(t) - Dic(t,-)] 

= A L x(t). 

Defining the velocity potential as before (Section 2) by 

<j)(t,x,z) = 0(0, x,z) + ( fi(a } x, z) dcr 

Jo 


in the notation 


and defining 


/i(V) 

x(t) = f 2 (t , •) 

/s(V) 

J^0(O ,x,z) = f 2 (0,x,z) 
-^-(/>(0,x,z) = f 3 (0,x,z) 


we see that (f)(t,x,z ) satisfies (2.1) and the boundary conditions. And 


0(0, x, z) 
0(0 ,x,z) 


can be specified arbitrarily within our far field conditions and boundary smoothness 
conditions and of course the differentiability conditions. The primary interest is 
however in the acceleration potential at z = 0+ definable through f 2 (t, x,z). 

Let t f)p(t,x,z ) denote the Possio solution for the velocity potential. This is a 
particular solution obtained by the “doublet-at-source” method. But this solution 
can be expressed in the form (5.3) for an appropriate initial condition. Let 

d<f>P 

dt 

Mt) = a -t 

d<j>p 
dz 


68 



Then we have: 


Xp(t ) = Dw(t,-) + S(t)(xp(0) — Dw(0, •)) — / S(t — a) Dw(cr, •) da. (5.6) 

Jo 

Let us define the Laplace transform: 

POO 

xp( A) = / e~ xt x P (t) dt, ReA>0. 

Jo 

Then we have 

x P ( A) = R(X,A)(x P {0) - Dw{ 0)) - R(X,A) (XDw(X) - Dw( 0)) + Du) (A) 

= Dw{ A) + R(A,Al)(xp(0) - XDw(X)) . (5.7) 

But since (see Section 3) 


<(>p(x,z, 0) = (j)p{x,z, 0) = 0 

we have that 

Xp(0) = 0 

and hence 

rr P (A) = (I — XR(X,A))Dib(X). (5.8) 

We have in (5.8) an alternate technique for finding Xp( A), alternate to the Possio 
equation. 

The acceleration potential defined by (2.22) can be expressed as a linear bounded 
transformation £: 

r/’(t) = Cx(t) 

and hence 

L’p(A) = 4> m { A) + £i?(A,/l)(ap(0) - XDw(X)) (5.9) 

which shows in particular that il’p(X) has an essential singularity along the entire 
negative real axis just as Dib(X) has but is otherwise analytic. 

Finally we note that 


£ (R(X, A)(x P {0) - XDw(X))) 

satisfies the Kutta-Joukowski conditions since both ipp(X) and do. But there 

are clearly many solutions of the field equations which do not. Indeed the question 
of characterizing those initial conditions which lead to solutions which satisfy the 
Kutta-Jouskowski conditions is open. 
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Finally we note that 


x(t) — x P (t ) = S(t)(x( 0) — x P ( 0)) 
or, 

x(t) = x P (t ) + S(t)(x( 0) — x P { 0)). 

Thus any two solutions differ by a term which depends only on the initial conditions 
assumed for the field. Since these conditions in a physical sense rarely are specifiable, 
the unsteady aerodynamics can never be determined completely. Note that 


\\S{t)x\\* 


so that there is no energy decay, within this theory! Of course solutions which do not 
satisfy the Kutta-Joukowski condition may be dismissed as “not physical.” 


General Representation for Velocity Potential 

The representation (5.1) can be recast as a representation for velocity potential. 
Thus: 


Theorem. 5.1 


Any solution f(t,x,z) of (2.1) with the flow-tangency and initial conditions can 
be represented as 

<p(t,x,z) = <p M (t,x,z) + 4> R (t,x,z) (5.10) 

where (p R (t,x,z) is the unique solution of (2.1) with a forcing term: viz.: 


d 2 4> R 

dt 2 




9 2 4>r 

dtdx 


<4(1 -M 2 ) 


d 2 (p R 

dx 2 


2 d 2 4> R _ d 2 (pu 
a °° dy' 2 - dt 2 


2 a oc M 


d 2 (p M 

dtdx 


(5.11) 


with homogeneous boundary conditions and initial conditions consistent with 
<p R (t,x,z) = <f)(t,x,z) - (t> M (t,x,z), t> 0. 

Moreover, (fi R (t,x,z) satisfies the Kutta-Joukowski conditions if (f)(t,x,z) does. 


Proof 


We calculate directly that 


£> 2 <Pr 
dt 2 


20,00 M 


d 2 f R 

dtdx 


d 2 f 

~dW 


+ 2 Cl'oo 1 1 


d 2 f 

dtdx 


d 2 4> M 

dt 2 


2a 0 0 M 


d 2 (f>M 

dtdx 
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Now 


d 2 4> 

dt 2 


2a 00 M 


d 2 (f) 

dtdx 


al(l~M 2 ) 


d 2 (f) 2 d 2 (j) 

dx 2 a °° dz 2 

d 2 (/)R 2 dVa 

dx 2 + a °° dz 2 


since 


(1 — M 2 ) 


d 2 4> M 

dx 2 


+ 


d 2 (p M 

dz 2 


0 


and (5.11) follows. Since both </>(•) and satisfy the flow tangency conditions, 

we see from (5.10) that 4>r(') satisfies homogeneous boundary conditions. The rest 
of the statement concerning 4>r(') is immediate. □ 


Remark 

Let (f>p(t,x,z ) denote the Possio solution. Then 

(j) R (t,x,z) = (/) P (t,x,z) - (p M (t,x,z) 
satisfies the Kutta-Joukowski condition and further 


0fl(O ,x,z) = 0 = 0 H (O ,x,z) 
if the downwash initial conditions are zero: 


iv a (0,x) = 0 = fr„((). ./•). \x\ < b 

and 

<t>R(t,x,z) 

is thus the response solely to the forcing function on the right in (5.11), with zero 
boundary as well as initial conditions. 
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